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The scattering of two hadronic objects at high energy is similar to a reaction-diffusion 
process described by the stochastic Fisher-Kolmogorov equation. This basic observation 
enables us to derive universal properties of the scattering amplitudes in a straightforward 
way, by borrowing some general results from statistical physics. 
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i/-) , 1. An event-by- event picture of high energy scattering 

O , 

We consider the scattering of a dipole of size r off a target |gg)-dipole of size r as a 
Ph! toy model for high energy scattering. In the QCD parton model and in the rest frame of 
o ,| the target, the projectile interacts through one of its quantum fluctuations \qqgggg- ■ ■) 
which are produced by QCD radiation. At high energy, these fluctuations are dominated 
by dense gluonic states, generically called "color glass condensate" [H]. 

It proves useful to view the partonic configurations as collections of color dipoles [ 
12] characterized by their transverse sizes r i,..., r n . Each of these dipoles may interact 
with the target according to the elementary amplitude T eZ (rj,r ) ~ a;f(r|/r>), where 
r< = min(rj,ro) and r> = max(rj,ro). The total amplitude for a given Fock state 
realization reads 

n 

T ( r , r o) = J2T d (ri,r ) . (1) 

i=l 

Roughly speaking, T is counting the number of dipoles within a bin of size 1 centered 
around p = log r 2 /r 2 , with a weight given by the interaction strength ot 2 s . 

Let us describe the rapidity evolution of a given Fock state. An increase of the rapidity 
of the projectile opens up the phase space for each dipole to split into two new dipoles 
(this is the dipole interpretation of gluon branching). The probability of such a splitting to 
occur is given by the Balitsky-Fadin-Kuraev-Lipatov (BFKL) kernel. The typical partonic 
configurations that drive the high energy evolution are those for which the newly created 
dipoles have sizes of the order of the size of their parent dipole. This means that the 
parton density grows essentially diffusively under rapidity evolution. 

The above picture is valid provided the interaction is not too strong, i.e. as long as 
T ^ 1. In the bins where the number of partons gets large and T reaches 1, Eq. (Q) and 
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the rapidity evolution itself are supplemented by nonlinear terms which tame the growth 
of T in such a way that the unitarity bound T < 1 be respected. At that point, the dipole 
picture itself breaks down. The exact mechanism of how unitarization is realized (which 
involves gluon recombinations and multiple scatterings) is still not fully understood, but 
such information is not required to get most of the properties of the amplitude, see 
Eqs. (1715)1 . 

The picture of high energy scattering that we have just outlined is essentially that of a 
reaction-diffusion process in a system made of N = I /a 2 particles [ |3] . 

So far, we have being focussing on one particular Fock state realization, which corre- 
sponds to one given event in an experiment. The amplitude T is the scattering amplitude 
for that given partonic state, which obviously is not an observable since partonic configu- 
rations are random and cannot be selected experimentally. The physical amplitude is the 
average of T over all partonic configurations accessible at rapidity Y: 

A(Y,r) = (T(r,r )) Y . (2) 

In the following sections, we will study the properties of A that can be deduced from our 
reaction-diffusion picture. 

2. Mean field approximation and the FKPP equation 

As a first step, we enforce a mean field approximation by neglecting the fluctuations of 
T (which are due to statistical fluctuations in the dipole number), in which case A = T. 
The rapidity evolution of A is well-understood in this approximation: it is given by 
the Balitsky-Kovchegov (BK) equation [ H] , which was recently shown to belong to the 
universality class of the Fisher-Kolmogorov-Petrovsky-Piscounov (FKPP) equation [|3J. 
The latter reads 

d- aY A = d 2 p A + A-A 2 , (3) 

where a = a s N c /n. The first two terms in the right handside stand for the growing dif- 
fusion of the partons, while the last nonlinear term tames this growth so that A complies 
with the unitarity limit. While the exact mapping between the BK and the FKPP equa- 
tions has not been exhibited, the universality class was unambiguously identified [|SJ. 
From the properties of the solutions of equations belonging to that class, it became clear 
that the two equations have the same large- F asymptotics, up to the replacement of the 
few parameters that characterize the diffusive growth, and that have to be taken from the 
BFKL kernel in the case of QCD. 

The FKPP equation admits asymptotic traveling waves (see Fig. [TJ, namely at large 
rapidity, A boils down to a function of a single variable 

A(Y,r)=A(r 2 Q 2 s (Y)) , (4) 

where the momentum Q s is called the saturation scale and is defined, for example, by 
the requirement that A(Y, 1/Q S (Y)) be some predefined number Aq ~ 0(1). The scaling 
law (j3J) is known as "geometric scaling" [E]. The rapidity dependence of Q s was found to 
be [00131 

\nrlQ 2 s {Y) = a^-Y , x ( 7 ) = 2^(1) - ^( 7 ) - 1>{1 - 7) , — = x'(7o) (5) 
7o 7o 
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where we have put back the parameters derived from the BFKL kernel x{l)- The last 
equation in (jSJ) defines 7o [Ej- Eq. © captures the leading- y behavior: two further sub- 
leading terms have been obtained recently [lEHS]- The asymptotic form of the amplitude 
is also known for r > 1/Q S (Y) [E1E], as well as the first corrections to the scaling (jlj) [ 

u. 

A very important property is that Q S {Y) is completely determined by the small- A tail 
which drives the evolution of the whole front. 
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Figure 1. Picture of the evolution of T{r) 
with rapidity (from left to right). Initially 
at Y = 0, T ~ T el (r,r ). When rapidity 
increases, more dipoles populate the typi- 
cal Fock state configuration, leading to a 
growing diffusion for T. Shortly after the 
unitarity limit T = 1 is reached, a traveling 
wave forms. 
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Figure 2. Different evolutions of a given 
initial condition over the rapidity ranges 
Y , 2Y , 3Y , 4Y (the bunches of curves 
gather 100 different realizations). One sees 
the dispersion of the front positions, which 
is of diffusive nature and thus of order 
A(lnr 2 Q^))~vW. 



3. Beyond the mean field: the "high energy QCD/sFKPP" correspondence 

The mean field description outlined above is justified when the number of dipoles of 
a given size is large. However, in regions in which the amplitude T is of order a 2 s , the 
statistical fluctuations dominate and the mean field approximation breaks down. Both 
the fact that the number of dipoles of a given size is finite, and that this number fluctuates 
from event to event resulting in fluctuations 5T ~ a s \/T, are neglected in the mean field 
description ©. Since the tail T C 1 determines Q s , we anticipate that these features 
have dramatic consequences. An evolution equation corrected for these effects may be 
written as 

d- aY T = d 2 p T + T-T 2 + a 8 yjT(l-T) V , (6) 

where rj is a Gaussian white noise. With respect to Eq. Q, we have added a noise term 
that accounts for the statistical fluctuations in the dipole number. Eq. © is the stochastic 
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Fisher-Kolmogorov (sFKPP) equation. 

Significant progress has been achieved in the last few years in understanding the prop- 
erties of its solutions [Ej. It has been realized that the discreteness of the number of 
dipoles induces large corrections to the saturation scale, namely [l9"| ITUl 13] 

«W = A-««y. (7) 

7o 2 In 1/aj 

On the other hand, the fluctuations in the dipole number result in a diffusive wandering 
of the saturation scale: different events have saturation scales which typically differ by 
(9(V aY). This is illustrated in Fig. El for a reaction-diffusion model that shares the 
essential properties of the full QCD problem (note that a s = 10~ 5 , therefore statistical 
fluctuations in the tail are not distinguishable on the plot). After having averaged over 
many events to get the physical amplitude (J2J), the wandering of the saturation scale 
causes a drastic breaking of geometric scaling. 

The variance of the wandering of the front has recently been quantified and this 
result enables us to derive the scaling law for A [ 3j that replaces geometric scaling (jlj) : 

A ( r,r)^(^mL). (8, 

Very recent work [El has aimed at matching these ideas with the so-called Balitsky- 
Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner ( JIMWLK) formulation of the 
color glass condensate [H], at the expense of a modification of the latter. We should stress 
that in any case, the results (J7I8|) directly stem from the physics of the parton model, 
and are presumably exact solutions of QCD: they are the leading terms in a large- y and 
small-a s expansion. 
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